function [g] = g_function_sim()
% creates difference of moments
%   inputs: 
%   outputs: g - difference between data and data_sim

% get simulated data
data_sim = simulate_data_answer();

% assign data mometns
u = 
tu = 
w = 
te = 
w2 = 

% solve for simulated data moments
u_sim = mean(data_sim(:,1));
tu_sim = 
w_sim = 
te_sim = 
w2_sim = 

% create g function should be (5 X 10000)
g = 

end

